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We theoretically and numerically analyze x-ray scattering from asymmetric-top molecules three- 
dimensionally aligned using elliptically polarized laser light. A rigid- rotor model is assumed. The 
principal axes of the polarizability tensor are assumed to coincide with the principal axes of the 
moment of inertia tensor. Several symmetries in the Hamiltonian are identified and exploited to 
enhance the efficiency of solving the time-dependent Schrodinger equation for each rotational state 
initially populated in a thermal ensemble. Using a phase-retrieval algorithm, the feasibility of struc- 
ture reconstruction from a quasi-adiabatically-aligned sample is illustrated for the organic molecule 
naphthalene. The spatial resolution achievable strongly depends on the laser parameters, the initial 
rotational temperature, and the x-ray pulse duration. We demonstrate that for a laser peak intensity 
of 5 TW/cm 2 , a laser pulse duration of 100 ps, a rotational temperature of 10 mK, and an x-ray 
pulse duration of 1 ps, the molecular structure may be probed at a resolution of 1 A. 

PACS numbers: 33.80.-b, 34.50.Rk, 61.Q5.cc, 42.30.Rx 



I. INTRODUCTION 

X-ray diffraction is a powerful method for investigat- 
ing structures of molecules. X-ray crystallography has 
become the standard tool for identifying the structure 
of large molecules and proteins [TH3J. Recent develop- 
ments in x-ray sources have opened new opportunities 
[4T113] for imaging membrane proteins and other macro- 
molecules that cannot be crystallized 14 . One possibil- 
ity is single-molecule imaging, where one molecule at a 
time is probed by an intense x-ray pulse [Trtll]. which 
subjects the molecule to severe damage [5] [TS]. A se- 
ries of diffraction patterns has to be collected and clas- 
sified according to the molecular orientation to get full 
structural information from randomly oriented molecules 
[9]. New iterative phase-retrieval algorithms for non- 
crystalline specimens have been developed to invert the 
diffraction data [SHU EEMEB] ■ 

Alternative approaches with laser-aligned molecules 
have been proposed, where an ensemble of molecules, 
rather than a single molecule, is exposed to x-ray pulses 
[U [5]. In this approach the radiation dose for each 
molecule is strongly reduced and stays well below the 
critical value for damage [J]. Three-dimensional infor- 
mation on the single-molecule structure can be gained 
from a well-aligned ensemble by capturing many two- 
dimensional diffraction patterns. The ability to accu- 
mulate the signal over a large number of x-ray pulses re- 
duces the radiation dose further. Notice that alignment, 
rather than orientation, of the molecules is expected to 
be sufficient for the reconstruction of the molecular struc- 
ture [iniOjl]. Nevertheless, various methods have been 
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proposed to orient molecules [501 [3T] by the linear Stark 
effect (22[, the AC Stark effect [23], or two-color fields 

Much theoretical and experimental research has been 
done on laser-induced alignment [35H7T]. With nonres- 
onant laser frequencies, an induced dipole moment can 
be created that couples back to the laser electric field 
and forces the molecule to be aligned. Earlier work has 
focused on one-dimensional alignment of the most polar- 
izable axis by using linearly polarized light [251456] . The 
demonstration of three-dimensional alignment [571463] . 
which requires an asymmetric-top molecule and an ellip- 
tically polarized pulse (or two linearly polarized pulses) 
opens the door for probing the three-dimensional struc- 
ture. The Coulomb explosion technique has been ex- 
ploited to detect three-dimensional alignment [SSI 
where rotational temperatures down to 1 K have been 
accomplished [72] . 

Depending on the duration of the laser pulse, T£, rel- 
ative to the rotational period of the molecule, r rot , the 
alignment dynamics can be classified into three distinct 
regimes. In the limits of adiabatic (t£ S> T rot ) and im- 
pulsive (jl <C T ro t) alignment, the dynamics can be de- 
scribed analytically [32J [33J . Impulsive alignment reveals 
the quantum mechanical nature of this process by show- 
ing alignment revivals after the laser pulse is over [34rt4T] . 
In the adiabatic limit, the alignment dynamics follow the 
laser pulse shape [55] , No analytic solution exists in the 
intermediate regime (tx ~ T ro t), and the time propaga- 
tion of the molecular ensemble in the presence of the laser 
pulse has to be performed numerically [33JEI]- Especially 
the quantum mechanical description of three-dimensional 
alignment [57] [5J2 [S3] has proved to be numerically ex- 
pensive [61] . 

A general theory of x-ray diffraction from laser-aligned 
symmetric-top molecules was developed in Ref. [52) . 
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Applications to adiabatically aligned molecules may be 
found in Refs. (52] [53]. Reference [40] describes cal- 
culations on x-ray scattering from impulsively aligned 
molecules, exploiting the alignment revivals for probing 
field-free molecular structure. It has been shown that, 
for a symmetric-top molecule containing a single heavy 
scatterer, a holographic algorithm can successfully recon- 
struct the molecular structure from an x-ray scattering 
pattern [53"] . 

In this article, we discuss x-ray scattering from an en- 
semble of rigid, asymmetric-top molecules aligned three- 
dimensionally by elliptically polarized light at finite ro- 
tational temperature. We point out the symmetries in 
the quantum mechanical theory of three-dimensionally 
aligned molecules that can be used to significantly re- 
duce the numerical time propagation. We restrict our 
analysis to electronic and vibrational ground-state con- 
figurations and neglect deformations. Our approach al- 
lows us to investigate x-ray diffraction from molecules 
in all three alignment regimes (adiabatic, impulsive, and 
intermediate). In order to probe molecular structure of 
gas-phase molecules by x-ray scattering, the degree of 
alignment must be rather high for sufficient resolution. 
This favors the adiabatic alignment approach |43] with 
low-temperature molecules [25j[55]. In Sec. [Tl| we present 
the theory and inherent symmetries of the Hamiltonian 
underlying x-ray diffraction from laser-aligned molecules. 
Section III focuses on the numerical implementation and 



computational efficiency, and basic ideas are presented 
of the phase-retrieval algorithm that is used for struc- 
ture reconstruction. In Sec. |IV[ the three-dimensional 
alignment dynamics and their impact on the reconstruc- 
tion are demonstrated using the example of the organic 
molecule naphthalene. We conclude with a discussion of 
the feasibility and limitations of x-ray diffraction from 
laser-aligned gas-phase molecules. Atomic units [73] are 
employed throughout, unless otherwise noted. 



II. THEORY 

The Hamiltonian for x-ray diffraction from laser- 
aligned molecules is [53] 



x-ray fields, and H^(t) and Hx describe the interactions 
of a molecule with the laser and x-ray field, respectively. 
The electronic and vibrational states of the molecule are 
omitted, since we assume that the molecule remains in its 
electronic and vibrational ground state throughout our 
discussion. The x-ray field and its interaction with the 
molecule is described in a quantized manner. The laser 
field is formulated as a classical field. The eigenstates of 
the noninteracting system, i.e., Hq = H YOt + Hxem, are 



with 



\JrM-{n}) = \JtM) ® \{n}) , 

Hxem \{rn}) = Efa \{m}) , 
H ro t | J tM) = E?*\JtM), 



(2) 

(3a) 
(3b) 



where \JtM) are the rotational eigenstates of an 
asymmetric-top molecule j57[ |?1] and \{n}) are the x- 
ray Fock states [52]. The density matrix of the whole 
system is 



Ptot 



x 



J2 J2 P{nt},{n2} W Jr 
JtM {m},{n 2 } 

X \^JrM;{ ni }(t)) (*JrM;{n 2 }(*)| 



(4) 



where the gas-phase molecules [75, 76i are described by a 
canonical ensemble [77M79] . wj T is the statistical weight, 
and i r i denotes the initial distribution of all the 



occupied field modes [SO] [81] . 

With including both interactions, H\,{t) and Hx, the 
states of Eq. |2]) are no longer eigenstates of the system. 
However, each state |\& j oTa M -{no}(^)) can be written as 

\Vj T M ;{no}(t)) = ( 5 ) 
JiT 1 M 1 .{n 1 } 

x^(*)]£;:£K I JmMr, {m}), 



where the expansion coefficients [C(t)] 
the initial condition 



Jo^M ;{n } StatlSt y 



Htot — H Iot + Hxem + Hj, (t) + H 



x- 



(1) 



\y\ ^JhoroMo-Ano} 



£joJi<Wi£MoMi^{no},{ni}-( 6 ) 



where H mt describes the field-free rotational motion of a 
molecule, Hxem is the Hamiltonian of the free quantized 



The equation of motion for the expansion coefficients 
in the interaction picture (subscript I) reads 



d 

"dt 



(7) 



x ^V^n +E? n2 } -E? ni} )t [c (t)] ^n Afx Am} 

*■ 1 y n J T Mo;{no} 

I 



In the following, we assume we found [CL,i{t)\ j^^S 



the solution for the laser-only problem, i.e., H Iot + H^(t). 



3 



Note that in the laser-only Hamiltonian no x-ray field is 
involved and we will drop the x-ray field indices in its 
solution. 



The interaction between the laser-aligned molecules 
and the x-ray field is taken into account by first-order 
perturbation theory. The solution of Eq. ([7]) becomes 



J J r M 



(8) 



J3T3M3 



J21-2M2 

The expectation values of interest can be calculated by 



( J 2 r 2 M 2 ; {n 4 }| ffx |JiTiM i; {n }) [C LJ (t')} 



JotoMo 



0(f) = Tr[6p tot (t) 



£ ^'},{n"} ( J 2T 2 M 2 ; {n 2 }| 6 I J m M i; {m}> 

{«2} : Wl}:{«'}^K'} 



(9) 



from which the alignment signals (cf. Sec. II D I and the 



x-ray scattering probability (cf. Sec. II E ) can be derived. 
We followed here the strategy that is laid out in Ref. 
[52] for x-ray diffraction from symmetric-top molecules 
and adapt it to the asymmetric-top case. The structures 
of Hxem arL d Hx are the same for symmetric-top and 
asymmetric-top molecules. The most dramatic changes 
for asymmetric-top molecules happen in the laser-only 
system. Therefore, we focus for the rest of our discussion 
just on the laser-only system. 

In the following subsections the structure of the field- 
free Hamiltonian H ro t an d the laser-molecule interaction 
matrix Hi,{t) will be investigated. For this purpose it 
is convenient to stay in the Schrodinger picture. In Sec. 
|II C| we return to the equation of motion of the laser- 
only system and point out the symmetries of its solutions 
[CL,i(t)] j^nf 1 ■ The measure of three-dimensional align- 



ment is described in Sec. |IID[ The theory section closes 
with exploiting the symmetries in the angular density 
distribution and the diffraction signal. 



A. Free asymmetric-top rotor 

Assuming that structural deformation of the molecule 
may be neglected, we treat the molecules as rigid rotors. 
The corresponding field-free Hamiltonian reads [71] 



= AJ 2 n +BJi 



A 



2 

A- 



B_ J2 



2C-A-B 



(10) 



B 



J+i + J-i 



where A, B, C are the rotational constants associated 
with the principal axes of inertia, J a , J c are the Carte- 
sian components of the angular momentum operator in 
the molecular frame; and J±i, Jo are the spherical ba- 
sis components. For symmetric-top (A = B 7^ C) and 
asymmetric-top (A ^ B =^ C, A ^ C) rotors, the angular 
momentum J and its projection on the space-fixed z-axis, 
M, are conserved. The angular-momentum projection on 
the molecular c-axis, K, is conserved only for symmetric- 
top molecules. A new quantum number r, replacing K, 
must be introduced for asymmetric-top rotors diagonal- 
izing the field- free Hamiltonian [cf. Eq. (3b)]. Note, the 



rotational energies E™* are independent of the quantum 
number M. Hence, each energy level E T f* is (2 J+ l)-fold 
degenerate. It is possible to express the asymmetric-top 
eigenstates | JtM) as a superposition of the analytically 
known symmetric-top eigenstates | J KM) [57], 



\JtM) 
(<j>,6, X \JKM) 



Y j a [ ^ ] T \JKM), 



(11) 



K 



2J+ 1 



D^i{^e,xl (12) 



8tt 2 



where the coefficients a$ r are re& l an( l 



9, x are the 

three Euler angles connecting the space-fixed laser frame 
(L) with the principal axes of inertia in the molecular 

reference frame (M). D^/L(ip,6,x) is the complex con- 
jugate of the Wigner D-matrix with angular momentum 
J. Unfortunately, it is not possible to find a general re- 
lation between the r and K classification [82"] . 

There are two symmetries incorporated between the 
symmetric-top (A') and asymmetric-top (r) classifica- 
tion. The first symmetry decouples states with even K 
from states with odd K, since H mt is a linear combination 
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of J"± x and Jg [Eq. (10)]. The state class with even K is 
labeled E, the one with odd K is labeled O. Asymmetric 
eigenstates inherit this separation and are superpositions 
of just even or odd K states. The invariance of H mt un- 
der the substitutions M — > —M and K — > —K leads to 
the second symmetry, the Wang symmetry [741177] . where 
asymmetric-top eigenstates decompose in symmetric and 
antisymmetric linear combinations, 

\JtM) = ^ a W[\JKM) + (-l) T \J(-K)M)],(13) 

K>0 

where (— 1) T denotes the parity of t, 
(-l) r = 



t symmetric in K, 
t antisymmetric in K. 



(14) 



Thus, overall there are four separate state classes, which 
are summarized in Table [i] 

TABLE I: Rotational classes of asymmetric-top rotors. 



E+ : 


\JtM) 


- E <4 J l 

K>0, K even 


E~ : 


\JtM) 


- E 4 J i 

K>0, K even 


0+ : 


\JtM) 


- E 4 J l 

K>0, K odd 


Q- : 


\JtM) 


= E 4 J l 

K >0, K odd 



] T [\JKM) + \J(~K)M)] 



With knowing the energy levels of iT ro t, it is possible 
to calculate the laser-only density matrix 



P m °\t) = E^I^mW) <<Wi)|, (15) 

JtM 



with the statistical weights, 



wj T 



9Jt 



e -E?;/kT 
Z(T) ' 



(16) 



where Z(T) is the partition function at temperature T 
and k is the Boltzmann constant. Every state is addi- 
tionally weighted by the nuclear spin statistical weight 
<?j T , which represents the number of allowed nuclear spin 
states for a given rotational state and can be derived 
from symmetry arguments that have to obey spin statis- 
tics [53]. The computation of g,j T for asymmetric-top 
molecules and in particular for naphthalene is discussed 
in Appendix [A] 



are the unit vectors of the major and minor polarization 
directions. The parameters e x and e 2 satisfy e x + t 2 z — 1 
and < e T < e,. The laser-molecule interaction reads 



H L (t) 



E 



pol j j 



(t) 



(18) 



2 



J 



2^ ^ 

J=0 M=-J 



(-1) 



where a po1 is the dipole-polarizability tensor and U(t) = 
E(t) £g) E(t) is the electric-field tensor. On the right- 
hand side of Eq. ( 18 1, the interaction is written first as a 



Cartesian tensor product and then as a spherical tensor 
product. Only spherical tensor components with J = 0, 2 
are non-zero and contribute to the tensor product, since 
both tensors are symmetric. (All components with J = 1 
are zero for symmetric tensors.) The J = component 
shifts all rotational energy levels by a state-independent 
amount and may therefore be dropped. 

The laser period, 2tt/uj, is typically several orders of 
magnitude smaller than the rotational time scale, r rot 
1/(A + B). Cycle averaging leads to a diagonal Uij, and 
only three spherical components remain nonzero, 



[U(t)] [ o ]L = ^(2-3e*)l(i), (19a) 



Ana 
27rae 2 x I(t). 



(19b) 



Hence, only the terms involving [aP 01 ]^ and [a po j+l x ' 
contribute to the laser-molecule interaction [Eq. (18)]. 
This holds for any molecule. However, the polarizabil- 
ity a po1 is a molecular property and is therefore most 
conveniently expressed in the molecular reference frame 
gner D-matrices provide the connection to 
the space-fixed components [59], 



\M 



K 



(0,0, x) KH 



poll [J]r 



K 



(20) 



When the principal axes of the polarizability tensor do 
not coincide with the principal axes of the moment-of- 
inertia tensor, the polarizability tensor is not diagonal 
in the molecular reference frame (cf. Sec. II A) and all 

[a pol ]^' A/ for J = 0, 2 may be nonzero. We will restrict 
our discussion to the case where both frames coincide. 
The remaining tensor components are [59) 



B. Laser— molecule interaction 

Three-dimensional alignment may be achieved by using 
an intense, nonresonant, elliptically polarized laser field, 

E(t) — y/%'KaI{t) (e x cos(wi) e x + e z sin(wt) e^j ,(17) 
where I(t) is the intensity of the laser field, a is the fine 
structure constant, u is the laser frequency, and e z ,e x 



v poll[2]j\,/ 

' Jo — 



v poll[2]j\f _ 
* J±2 — 



2a po1 



"6h u aa 



yPol 



V6 

pol 



(21a) 
(21b) 



Using Eqs. (19) - (21), the matrix elements of the 
laser-molecule interaction operator with respect to the 
symmetric-top eigenstates read 



5 



(JKM\H L (t) \J'K'M' X 



1 2J + 1 



2 V 2 J' + 1 

x ( [a pol } l a Ul ( J, K; 2, 0| f, K 1 ) + [a^]f M \(J,K; 2, 2| J'K') + (J,K; 2, -2| J', K') 



(22) 



x [ |7 C*)]o 2] ^ ( -V/ : 2. ()| ./'. J /'; - [f ( / )] " ( J, M; 2, -2| J', M') + (J, M; 2, 2| J', M') 

I 



where the matrix elements of the Wigner D-matrices have 
been expressed in terms of Clebsch-Gordan coefhcients 

pa cm. 

Three-dimensional alignment of asymmetric-top 
molecules conserves neither K nor M . The asymmetric- 
top rotor breaks the x-symmetry (K conservation), and 
elliptically polarized light breaks the 0-symmetry (M 
conservation). However, there are remaining symmetries 
in the laser interaction that can be employed, 



(JKM\H L (t) \J'K'M') 
= 0, K' - K, M' - M 
■ J- J' 



{±2,0}, 



(23a) 



(-1) J - J (JK{-M)\H L {t)\J'K\-M')), (23b) 



= (-1)' 7 - 7 ' (J(-K)M\ H L {t) \J'{-K')M') . (23c) 

This follows from properties of the C lebsch-Gordan coef- 
ficients [71] . As a result of Eq. ( 23a I , states with even K 



(M) remain separate from states with odd K (M). The 
consequences of Eqs. (23b) and (23c I will be discussed 
in the following sections. 



but breaks the Wang symmetry [cf. Eq. ( 13 )] . The 
indices from the asymmetric-top basis transform to the 
symmetric-top basis according to 



a^Url[C L At)\V T u -(26) 



J't'M' 



The symmetric-top representation is especially favorable 
for the laser-molecules interaction; the field-free propa- 
gation is naturally expressed in the asymmetric-top basis. 



D. Measure of alignment 

The degree of three-dimensional alignment can be 
characterized in terms of the quantities cos 2 #z m , where 
9i m is the angle between the space-fixed axis / and the 
body-fixed axis m. The following relations hold among 
the cos 2 6im [551 : 



cos2 °lm 



^cos 2 6» iom = l Vm , l , (27) 



C. Equation of motion 

In the laser-only system, the equation of motion [Eq. 
([7])] for the initially populated state | J KM) reduces to 



4 [ClA^jkm = E < JVM 'I I J ^ M ^ ( 24 ) 

' JiTiMt 



xe^V-^W* [C L j(t)]j\ 



JKM ' 



where the interaction matrix elements are expressed in 
the symmetric-top basis and the symmetries from Eqs. 
(23) can be used. 



The symmetries of the rotational eigenstates | JtM) 
[cf. Eq. (13) and Table [i] in combination with the sym- 
metries of Hl {t) [cf. Eqs. (23)] pass on to [Cl.j(£)] jkm 
such that 



= (-iy- J [c L , / (i)]f ( ( :^'if'(25b) 

The laser-molecule interaction preserves the separa- 



tion of even or odd K and M states [cf. Eq. (23a)] 



where five of the six relations are independent. These 
relations reduce the number of independent cos 2 0i m to 
4. The matrix elements (JKM\ cos 2 6 lm \J'K'M') of one 
set of independent cos 2 8i m (I £ {x 7 z},m e {a, c}) are 
given in Appendix |b| The symmetries of cos 2 9i m are the 
same as for H^t), 



(J KM | cos 2 e hn \J'K'M') 
= 0, K' - K, M' - M (£ {±2, 0}, 



(28a) 



= (-1 
= (-1 



j-j> 
j-j' 



JK(-M)\cos 2 6i m \J'K'(-M')) , (28b) 
(J(-K)M\ cos 2 9 lm \J'(-K')M') . (28c) 



The restriction imposed on K' — K by Eq. ( 28a ) makes 
it attractive to store the matrix elements of cos 2 Q\ m in 
the symmetric-top basis; the same is true for Hz,(t). 

The ensemble-averaged expectation values at time t 
are 



(cos 2 e lm ) (t) 

(cos 2 B lm ) JrM {t) 



JtM 



w jT (cos 2 lm ) jTM (t), (29a) 



E E 

J1K1M1 J 2 K 2 M 2 

x UiKxMx] cos 2 9 lm \J 2 K 2 M, 



(29b) 



x[C*lM]£m 



6 



with 



VMjrU = 2^ a T>K> e J > T \ C L,l{t)] JtM (30) 



Making use of the symmetries in Eqs. ( 25 ) and ( 28 ), each 
(cos 2 Qim) M (t) fulfills the relations 



(cos 2 e hn ) JKM (t) = (cos 2 e lm ) 

(co S 2 e lm ) jTM (t) = (cos 2 6 lm ) (t). (31b) 



Eq. (31a I cannot be used for asymmetric-top molecules, 



since K is not a good quantum number. However, this 
equation can be exploit for symmetric-top and linear ro- 
tor calculations. 



E. Angular probability distribution 

For the calculation of the x-ray scattering probability 
t£ the solution for the full Hamiltonian [Eq. nib] has to 
be known. The laser-only problem is solved as aescribed 
in Sec. II C| The x-ray interaction is treated in first-order 
perturbation theory as indicated in Eq. (J8|. Further- 
more, the following assumptions of the x-ray pulse are 
made: (1) The coherence time of the x-ray pulse is signifi- 
cantly larger as the rotational time scale of the molecules. 
(2) The bandwidth of the x-ray pulse is much larger than 
any rotational transition energy. (3) The angular spread 
of the x-ray pulse is considerably small. In Ref. [52], the 
detailed derivation of ^ is given with the final result 



dP dcr t h 



5(0), 



(32) 



where dat^/dfl is the Thomson scattering cross section, 

S(Q) = J dt j x {t) JJJ d(f>d8d X sin 6 (33) 
Xpfo,0, X ;t) |F mol (Q, X ) | 2 

is the diffraction signal, and jx{t) is the x-ray flux. The 
molecular information in the diffraction signal is con- 
tained in the molecular form factor, 



F mo} (Q,<t>,e, x ) 



d 3 r M p(r M )e- lQ [ R ^ e ^ r ^4) 



which is a function of the molecular orientation, i.e., 
Euler angles, and the momentum transfer Q given in 
the space-fixed frame (L). The integration d 3 rM is done 
in the molecular rest frame (M) with p{tm) being the 
electron density of the molecule. The rotation matrix 
R((j),6,x) transforms the vector tm into the space-fixed 
frame. The angular probability distribution, 



w jT \(4>,o,x\* Jtm(*))| ■ 

JtM 



(35) 



is linked to the Wigner D- matrices k ( 
Eqs. (fill and p|. 



0, x) through 



From the fact that the molecules are aligned rather 
than oriented, symmetries additional to the ones of 
d m,k(^ ^X) EU enter into p(4>,9,x;t). These addi- 
tional symmetries originate from Eqs. ( |25[ ) and the sep- 
aration of the rotational classes E ± and ± throughout 
the alignment process (cf. Sec. II B). In terms of Euler 



angels, the symmetries of p((f>, 9, x; t) are: 



p(0,0,x;*) 



p{<t> + A", 6, x; t) = p(-<f>, ir-9, X + ir; t), (36a) 
p{4>, 6, X + TT;t) = p{^ + ir,n-e, -x; t). (36b) 



The first two symmetries [Eq. (36a)] correspond to C2 
rotations about the space-fixed axes z and x, respec- 
tively. The last two symmetries [Eq. (36b)] correspond 
to C2 rotations about the body-fixed axes c and a, re- 
spectively. In Eq. (34), these symmetries translate 
into a replacement of R(4>, 0, x) by C2 tX / z R{<j>^-,X) an d 
R(4>, 0, x)^2,a/ct respectively. By letting the C 2 , x /z rota- 
tions act on Q rather than on t*m, the S(Q) symmetries 
are found, 



S(~Q x ,Q y ,Q z ), (37a) 
S{Q x ,-Q yi Q z ), (37b) 
S(Q x ,Q y ,-Qz)- (37c) 



Friedel law [55], i.e., F*(Q) = 
S(—Q), has been used. The sym- 



Additionally, the 
F(-Q) S(Q) = 
metries of Eq. ( 36b ) can be used to reduce the integra- 
tion range of the Euler angles in Eq. ( 33 ) but only the 
symmetries of Eq 



(36a 



survive the integration, which 
are expressed in Cartesian coordinates in Eqs. (37). Re- 
gardless of their internal structure, Eqs. (37) hold for 



all aligned molecules. Note that alignment does not dis- 
tinguish between parallel and anti-parallel orientations. 
Thus, even for perfect alignment, i.e., (j) = 6 = X = 0, 
four distinct molecular orientations contribute incoher- 
ently to the diffraction signal S(Q) [19 . 



III. NUMERICAL METHODS 

Earlier work has addressed the problem of numeri- 
cal efficiency in the computational treatment of three- 
dimensional alignment of asymmetric-top molecules (61j . 
In Sec. |III A[ we describe numerical techniques and sym- 
metry arguments we have implemented to decrease the 
numerical effort. As far as we are aware, these specific 
points have not been discussed earlier in the literature. 
In Sec. |IIIB[ we explain the phase-retrieval method we 
employ to reconstruct molecular structures from x-ray 
scattering patterns. 
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A. Alignment dynamics 

The density matrix of a canonical ensemble is a sum 
of density matrices | JtM) wj t ( JtM\ [cf. Eq. (15)]. 
Analogously, all observables can be written as a sum 
of independent contributions, which may be calcul ated 
separately. This was shown for (cos 2 6>/ TO ) [Eq. (29)] 
and p(cj),6,x]t) [Eq. (35)]. Each contribution has a 
well-defined behavior under the substitution M — > —M, 
which can be used to avoid the calculation of \^j T M(t)) 
for M < 0. The same is true for the quantum number 
K, but unfortunately K is not a good quantum num- 
ber for asymmetric-top molecules, and the wave function 
\^jTAi(t)) has to be known to calculate the expectation 
values. However, \^fj T M(t)) can be written as a super- 
position of \^jKM(t)) [cf. Eq. ([ll])], and by using the 
same argument as for M, only jWjiof (t)) for K > 
are necessary to build all \"9j rM (t)) [cf. Eq.(13)]. This 



makes it attractive to propagate \^jKAi(t)} rather than 
\^ jTAi(t)) ■ Taking both symmetries together only the 
states \^jKM(t)) for < K, M < J have to be propa- 
gated to understand the full system response, which can 
save up to a factor 4 in computational effort. 

The decoupling between even and odd K/ r-states and 
M-states also enhances efficiency and does not get de- 
stroyed by the presence of the interaction Hl (t) . There- 
fore, many [Cl^Yjkm remain zero throughout the 
alignment process. This holds for symmetric-top as good 
as for asymmetric-top rotors, since each \Jt M) c an be 



classified into the class E ± or ± (cf. Sec. II A). The 
benefit is not just a speed-up by a factor 4; also the mem- 
ory requirement to store [Ch{t)\jKM > the matrix i?z,(t), 
and the matrices cos 2 6*/ m is reduced by a factor 4. Mem- 
ory size can become an issue on PCs when high J states 
are not negligible. 

A significant time factor in the calculation, besides 
the time propagation, is the computation of p(cj),6,x',t). 
By employing all symmetries described in Eq. ( 36 1 , the 



range of Euler angles with nonredundant information is 
reduced by a factor of 16 compared to the entire domain 
of Euler angles. 

Besides the physical symmetries that can be retrieved 
from the Hamiltonian, there are two numerical aspects 
that may improve the propagation speed. First, the 
numerical time propagation is commonly done by the 
fourth-order Runge-Kutta method, where it is impor- 
tant that the change in the wave function per propa- 
gation step is small and stays in the convergent region 
[86] . Yet too small step sizes quickly lead to numerical 

inefficiency. For our problem the change in [CL(f)\ J jKM 
is directly proportional to the product I(t)dt, where dt is 
the propagation time step. With decreasing I(t), dt can 
be chosen larger without leaving the convergent region. 
Variable step sizes are, therefore, important and improve 
the propagation efficiency further. Numerical tests have 
shown that the optimized program runs faster by a fac- 
tor between 3 and 4, where we assumed a Gaussian laser 



pulse (•entered at to with a full width at half maximum of 
tl- Outside of the laser pulse, i.e., t ^ [to — 3tl , to + 3tl] , 
the molecules are treated as field-free, and the propaga- 
tion is performed analytically. 

The second improvement takes place at the equation of 
motion [Eq. |24|]. To solve the equation numerically, it 
has to be discretized in time. We chose the fourth-order 
Runge-Kutta method as our discretization method. Only 
the contribution of H^t) is approximated when the dis- 
cretization in time in the equation of motion is made in 
the interaction picture. The time propagation of H mt is 
analytically exact such that in the field-free case the nu- 
merically calculated solution matches the analytical re- 
sult. To improve the efficiency in solving the equation 
of motion, we reduce the need for repeated calculating 
of e n j ; ' , which goes with N for the propagation 
of all [CL,i(t)] jk m , where N is the number of rota- 
tional states involved in the propagation. The evaluation 
of exponential functions does not occur when the equa- 
tion of motion is discretized in the Schrodinger picture 
[27l l28l [57] . As a consequence, the propagation of the 
field-free part is now discretized, which has two major 
limitations: (1) the field-free propagation is not exact 
even in the field-free case, where the analytic solution is 
known, and (2) the propagation step size depends, in ad- 
dition to the laser-molecule interaction strength, on the 
highest rotational energy E T f*. In order to use the ad- 
vantages of both pictures, we discretize the equation of 
motion in the interaction picture and transform it after- 
wards back into the Schrodinger picture by reformulating 

the propagation in terms of [Cl^Yjkm • The nna l dis~ 
cretized equation of motion is: 
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where the indices k are shortcuts for the sets 
of asymmetric-top quantum numbers, [-ffi(t)]^ — 
(j\H L (t)\i), e -^ E T dt l 2 = e" 1 ^, and intermediate so- 
lutions are 
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Since e~ tlf,j depends only on dt and not on the time t 
itself, it needs to be evaluated only once at the beginning 
of the propagation. During the time propagation, the 
evaluation of the exponential function is not necessary as 
long as dt does not change. In our simulations, we gained 
a factor 2 to 5 depending on the simulation parameters, 
when we employed Eqs. ( 38 1-( 40 ) . Note, in the field-free 



limit the time propagation in Eq. ( 38 ) is analytically 
exact. 

Overall, the use of the symmetries and numerical tech- 
niques described allows us to simulate three-dimensional 
alignment of asymmetric-top molecules from the impul- 
sive to the adiabatic regime up to two orders of mag- 
nitude faster. In the case of quasi-adiabatic alignment 
of naphthalene at 1 K, the total computation time was 
2457 h (41 days) with a PC (CPU: 3 GHz). The calcu- 
lation of S(Q) took approximately 8% of the total time 
(0, 9, x, i-grid points: 25 x 30 x 30 x 11). Without exploit- 
ing the symmetries for the Euler angles, this calculation 
would be almost 1.5 times longer than the time propa- 
gation itself. A detailed list of the physical parameters 
used for the computations is given in Sec. |TV} 



B. Phase-retrieval algorithm 

For imperfectly aligned molecules, the diffraction pat- 
tern is an incoherent average of single-molecule diffrac- 
tion patterns of different orientations. In the limit that 
a high degree of molecular alignment is attained, the ob- 
tained pattern can be approximated as a single-molecule 
coherent diffraction pattern. (This is true for naphtha- 
lene, since orientation and alignment are equivalent for 
this molecule.) One may thus retrieve structural informa- 
tion from the single-molecule electron density map, which 
is the Fourier transform of the single-molecule scattering 
form factor, F(Q). Since the diffraction pattern provides 
only |F(Q)|, we need to obtain the associated phase be- 
fore we can recover the molecular structure. There are it- 
erative numerical algorithms that permit reconstructing 
the phase directly from the intensity data [T7J [T51 1571 - 
1101 j . These algorithms require intensity data sampled at 
twice the Nyquist frequency. Successful structural recon- 
struction using these algorithms has been demonstrated 
with experimental data [7| P^ l HTMTTU] . 

Here we use the hybrid-input-output (HIO) algorithm 
[8"9Tl9"T] . which involves iterative Fourier transformation 
back and forth between the object and Fourier do- 
mains. A solution is found when the known constraints 
are satisfied in both domains. We begin by obtaining 
an initial estimate of the object electronic density via 
an inverse Fourier transformation of the form factor, 
F(Q) = |F(Q)|e l< ^, which is obtained by assigning a 
random phase, </>q, to the measured modulus, |F(Q)|. 
The random phase 4>q is chosen such that Friedel's law, 
F*(Q) = F(—Q), is satisfied. With this estimate, we ini- 
tiate an iterative four-step algorithm, in which the k-th 
iteration is given as follows: 



1. Fourier transform of the object electron density, Pfc(r), 
to obtain Fk(Q). 

2. A Fourier domain operation (FDO) is applied to 
Fk(Q) to obtain F' k (Q) that satisfies the Fourier con- 
straint. In our FDO, the modulus of F' k (Q) is set to be 
the measured modulus, |F(Q)|, and the phase of F k (Q) 
is the phase of Fi c (Q). 

3. Inverse Fourier transform of F' k (Q) to give p' k (r). 

4. An object domain operation (ODO) is applied to 
p' k (r) to get a new estimate of the object electron density, 
Pk+i{ r ), that satisfies the object constraint. Our ODO 
is given as 




, for r S S and p' k (r) > 
, otherwise, 

(41) 



where /3 is chosen to be 0.9 [93l EHJ [100], and S is a 
pre-defined support of the object. 

In order to obtain a correctly reconstructed object, a 
support S of good quality is needed [jJTJ. In our algo- 
rithm, the support iS is changed dynamically throughout 
the HIO algorithm via the Shrink-wrap (SW) procedure 
[93] . The inverse Fourier transform of the scattering in- 
tensity, \F(Q)\ 2 , equals the auto-correlation function of 
p(r). Treating the auto-correlation function as a distri- 
bution function, the initial S is chosen as a region cen- 
tered at the mean of the auto-correlation function with a 
spatial extension of two standard deviations. This choice 
of iS is fixed during the first 2000 iterations before apply- 
ing the SW procedure periodically after every 200 iter- 
ations of the HIO algorithm to obtain a new S. In the 
SW procedure, the modulus of the object, |pfc(r)|, is con- 
volved with a Gaussian of width a . The new S is then 
selected as the region for which the value of the convolved 
function is above a threshold of 20% of its maximum [931 . 
The initial width of the Gaussian is chosen to be 2.5 A 
and is shrunk linearly to a minimum of 0.25 A after 400 
iterations of the SW procedure. Using the last updated 
S, an additional 200 iterations of the HIO algorithm are 
performed. 

Figure [T] illustrates the application of the phase- 
retrieval algorithm described to an x-ray scattering 
pattern calculated for perfectly aligned naphthalene 
molecules with one set of initial random phases. X-ray 
scattering patterns in this article are based on the as- 
sumption that the molecular electron density equals the 
sum over spherically-averaged atomic electron densities 
|85j . Only the carbon atoms are visible in the recon- 
struction in Fig. |ljb) , since x-ray scattering from carbon 
is much stronger than x-ray scattering from hydrogen. 

Note that the real-space pixel size in each dimension is 
given by 7r/Q max , where Q max is the maximum momen- 
tum transfer for which diffraction data are available. In 
the case of perfect alignment, the real-space resolution, 
d vcs , attained from the reconstruction is two times the 
real-space pixel size. (In Fig. [I] d rcs =0.25 A.) However, 
in the case of imperfect molecular alignment, g? ics de- 
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FIG. 1: (Color online) (a) The diffraction signal S(Q y = 0) 
for naphthalene perfectly aligned in the xz plane. The spatial 
resolution is 0.25 A. 5(0) = 68 2 is the number of electrons 
squared, (b) The reconstructed structure. 



pends on Q max and Q C oh , where Q co h is the range of use- 
ful diffraction data for which the assumption of coherent 
scattering holds and the effect of incoherent averaging is 
small. In fact, we find that including diffraction data be- 
yond Q C oh can diminish the quality of the reconstructed 
object. 



IV. RESULTS 

In this section, we present numerical results based on 
the theory and the numerical strategies summarized in 
Sees. [II] and [Hi] We demonstrate quasi-adiabatic, three- 
dimensional alignment of the organic molecule naphtha- 
lene (CioHg). The experimental parameters that we are 
using are based on previous work [56]. X-ray energies 
and fluxes that are used in our discussion are accessi- 
ble at the Advanced Photon Source at Argonne National 



Laboratory. 

After investigating the alignment dynamics, we focus 
on the diffraction signal and structure reconstruction. 
The rotational temperature and the x-ray pulse dura- 
tion impact the effective alignment and limit the struc- 
tural information stored in the x-ray scattering pattern. 
A phase-retrieval algorithm is used to reconstruct the 
structure. In this context, the degree of alignment and 
its impact on the spatial resolution are discussed. 



A. Three-dimensional alignment 

Naphthalene is a planar molecule. The molecular ref- 
erence frame is chosen such that all atoms lie in the 
ac plane and the b axis is perpendicular to it. The 
choice of axes as well as the structure of naphthalene 
are shown in Fig. [2] The rotational constants of naph- 
thalene are A = 0.041 cm -1 , B = 0.029 cm -1 , and 
C = 0.104 cm -1 [m][l~l"2]. The polarizability constants 
vP Ql = 121.4 0%, alf = 63.2 af ™ A ~ po1 



are a* 



D, and o£f = 163.9 ag 



|113j . The nuclear spin statistical weights gj T for naph- 
thalene are given in Table [TT} (We assume that the carbon 
nuclei are 12 C isotopes, and the hydrogen nuclei are pro- 
tons.) The derivation of gj T is outlined in Appendix [A] 




FIG. 2: (Color online) The structure of naphthalene (Cio-ffs)- 
All atoms lie in the ac plane of the body-fixed reference frame. 
The coordinate system is shown in the lower left corner. The 
cones symbolize the resolution limit of the carbon atoms due 
to imperfect alignment. 



[121] In contrast to spectroscopic convention, the rotational constants 
A, B and C are not ordered according to magnitude. We re- 
order the rotational constants such that for perfect alignment 
the molecular frame (a, b, c) coincides with the space-fixed frame 
(x,y, z). 
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TABLE II: The nuclear spin statistical weights gj T of naph- 
thalene classified by the rotational symmetry classes. 



Rotation class 


E+ 


E~ 


o+ 




gj T (J even) 


76 


60 


60 


60 


g jT (J odd) 


60 


76 


60 


60 



Three-dimensional alignment of naphthalene is ob- 
tained by employing elliptically polarized light. The z 
and x axes of the space-fixed frame are defined by the 
major and minor axes, respectively, of the ellipse that 
characterizes the elliptically polarized light. Following 
Ref. [53], we choose the ratio between the field compo- 
nents using 



Q' 



pol 



vPOl 



a. 



pol 



0.422, 



(42) 



'bb 



which maximizes the three-dimensional alignment of 
naphthalene. Perfect three-dimensional alignment with 
respect to the space-fixed frame is achieved, when 
(cos 2 8 x a ) — (cos 2 9 y b) = (cos 2 ZC ) = 1< Random orien- 
tation corresponds to (cos 2 0i m ) = 1/3 for all angles. 

The typical distance between two neighboring carbon 
atoms in naphthalene is 1.4 A. In order to resolve the 
atomic structure, the resolution must be smaller than 
this value. Imperfect alignment limits the resolution. To 
build a connection between the resolution and the align- 
ment in terms of (cos 2 0j m ), it is helpful to view the angle 
0i m as an opening angle within which the residual motion 
of the atoms takes place (cf. Fig. [2]). Hence, the length 
characterizing the smallest resolvable structure is 



2Ry/l - (cos 2 01, 



(43) 



where R is the linear dimension of the molecule (mea- 
sured from the center of mass of the molecule). The 
degree of alignment needed for a resolution of 1 A for the 
outermost carbon atom (R = 2.5 A) is (cos 2 9i m ) = 0.96. 
For a given resolution, smaller opening angles 9i m are 
required when the molecules become larger. In the adi- 
abatic regime at sufficiently low temperature the maxi- 
mum alignment for a linear rotor (A = B , C — 0) is given 
by (cos 2 6 ZC ) = 1-^/45/7,7= VMnal [a? 01 ] {f 1 [3H33J. 

The polarizability [a pol ]{) 2 ^ is approximately proportional 
to R 3 , and B scales approximately as R~ 5 [THIEO!- The 
smallest resolvable dimension d co h is thus proportional 
to 1/R and might be expected to decrease with increas- 
ing molecular size. Note, however, that the temperature 
required to suppress thermal effects (kT/B <C 1) also 
decreases as a function of the molecular size. Addition- 
ally, in order to remain in the adiabatic regime, the laser 
pulse duration must increase with increasing molecular 
size. In practice this means the laser intensity will de- 
crease. For large molecules it is more realistic to con- 
sider the high temperature limit (kT/B 1), where 
the degree of alignment is a competitive interplay be- 
tween rotational temperature and coupling strength 7; 



more precisely, (cos 2 6 ZC ) = 1 — y/nkT/j [25] 133] ■ By us- 
ing the same scaling arguments, we find c? co h oc y/kTR. 
The expected resolution is now reversed and increases 
with molecular size and rotational temperature. Conse- 
quently, larger molecules can be resolved less precisely. 
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FIG. 3: (Color online) Alignment dynamics of three axes of 
naphthalene at rotational temperatures T = 10 mK (solid) 
and T = 1 K (dotted). The pump laser shown in gray has a 
peak intensity of / = 5 TW/cm 2 , a width of tl = 100 ps, and 
an ellipticity of e 2 /e 2 = 0.422. The upper panel pictures the 
alignment dynamics on the time scale of the pulse duration. 
The lower panel is a close-up view of the region highlighted 
in the upper graph. At t = 20.75 ps, nonadiabatic oscillation 
enhances the alignment for all axes at T = 10 mK and is ideal 
for a 1-ps x-ray probe pulse. 



In Fig. [3j the alignment dynamics of naphthalene 
at 10 mK and 1 K, respectively, are shown. We as- 
sumed a Gaussian laser pulse with a peak intensity of 
7 = 5 TW/cm 2 and a pulse duration of 77, = 100 ps 
(FWHM). At a rotational temperature of T = 10 mK, 
almost all naphthalene molecules are in the rotational 
ground state. The rotational period of naphthalene at 
such a low temperature, r rot sa 1/(A + B) ps 476 ps, is 
of the order of tl = 100 ps, suggesting that the align- 
ment dynamics are quasi-adiabatic [55) . This is consis- 
tent with the observation that the overall alignment fol- 
lows the laser pulse shape and, in addition, clear nona- 
diabatic features (fast oscillations) are visible. One can 
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make use of the nonadiabatic behavior by probing the 
molecules with an x-ray pulse of 1 ps duration, which is 
fast enough to resolve the nonadiabatic oscillations. For 
instance, at t — 20.75 ps the alignment is transiently en- 
hanced in all three dimensions (cf. Fig. [3]). At T = 1 K, 
the rotational period of naphthalene is small in compar- 
ison to T£, so the alignment dynamics are adiabatic and 
fast oscillations are significantly suppressed. As a con- 
sequence of the increased thermal motion, the maximum 
degree of alignment at T = 1 K is clearly reduced. 

The diffraction signal S(Q) [cf. Eq. ((33))] collected 
over the x-ray pulse duration reflects a pulse-averaged, 
effective alignment 



off 



dtj x (t) (cos 2 9 lm ){t), (44) 



where jx{t) is the normalized x-ray flux. We assume a 
Gaussian temporal envelope for the x-raypulse, with a 
full width at half maximum of tx- In Fig. Pfl (cos 2 9im) eS 
is shown as a function of Tx- We may conclude from 
Fig. [4] that the effective alignment of the body-fixed axes 
decreases rapidly when the x-ray pulse duration is longer 
than the laser pulse duration. An enhancement in the 
effective alignment is visible for Tx ~ 1 ps at T — 10 mK, 
where nonadiabatic oscillations are not suppressed and 
can be resolved by the x-ray pulse. 
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FIG. 4: (Color online) The effective alignment of the body- 
fixed axes is shown as a function of the x-ray pulse duration, 
T Xy for T = 10 mK (solid) and T = 1 K (dotted). The 
pump laser has a peak intensity of I = 5 TW/cm 2 , a width 



of tl — 100 ps, and an ellipticity of e 2 /e 2 
pulses are centered at t — ps. 



0.422. The x-ray 



laser field is applied in this direction. The strong align- 
ment of the b axis translates to well-aligned naphtha- 
lene molecules in the polarization plane of the laser (xz 
plane). Within the plane, the body-fixed c axis is always 
more strongly aligned in the z direction than the a axis 
is aligned in the x direction. (Recall that a?" 1 > ag" 1 
and e z > e x .) For higher temperatures, the body-fixed c 
axis is the most strongly aligned axis, since the alignment 
of the b axis is affected by the more rapidly decreasing 
alignment of the a axis. 



1 



0.96 ■ 



% 

A 
E 
a 



0.92 



0.88 





1 




cos 2 ^ 






— 






cos 2 e IH — 












t x -1 00 ps 


■ 





0.25 0.5 

temperature [K] 



0,75 



FIG. 5: (Color online) The effective alignment of the body- 
fixed axes is shown as a function of the rotational temperature 
for tx — 1 ps (solid) and tx ~ 100 ps (dotted). The pump 
laser has a peak intensity of / = 5 TW/cm 2 , a width of tl = 
100 ps, and an ellipticity of e 2 /e 2 = 0.422. The x-ray pulses 
are centered at t = ps. 



As may be seen in Fig. (4) the effect ive alignment 
(cos 2 0) s of all molecular axes is at least 0.96 for 10 mK 
naphthalene molecules probed by x-ray pulses shorter 
than ~ 20 ps. Hence, the smallest resolvable dimension 
limited by residual pendular motion is d co h — 1 A. Only 
the rotational ground state is initially occupied, and the 
minimum intensity required to accomplish 1 A resolution 
is I = 5 TW/cm 2 . Dissociation and ionization play only 
a minor role at this intensity, but become important at 
higher intensities |114j . With an effective alignment of 
(cos 2 #xa) eff = 0.88 at 1 K and tx = 100 ps, the smallest 
resolvable dimension G? co h is 1.73 A, which is larger than 
the distance between neighboring carbon atoms. 



When tx < tl, the influence of tx is rather weak 
and the degree of alignment is mainly affected by the 
rotational temperature, T. Figure [5] shows the effective 
alignment of the molecular axes as a function of T for 
tx = 1 ps and tx — 100 ps. Below 0.25 K, the best 
aligned axis is the body-fixed b axis in the space-fixed 
y direction for tx < 100 ps, despite the fact that no 



B. X-ray diffraction patterns 

In Fig. (6) diffraction patterns S(Q) of three- 
dimensionally aligned naphthalene are shown for perfect 
alignment as well as for imperfect alignment at differ- 
ent temperatures and x-ray pulse widths. The laser pa- 
rameters are the same as for Fig. |3] The laser polar- 
ization plane is assumed to be the xz plane so that the 
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molecules are aligned as illustrated in Fig. [2j The diffrac- 
tion patterns in Fig. [6] are two-dimensional planar slices 
through the three-dimensional Q space, taken at Q y = 0. 
(Due to the curvature of the Ewald sphere [55], experi- 
mental scattering patterns do not correspond to exactly 
planar slices through Q space, but this is of no conse- 
quence here.) The signal strength falls rapidly for high 
Q. Therefore, to highlight the structure at high momen- 
tum transfers, the scattering intensities are shown on a 
logarithmic scale. The effective alignment in Fig. [6] de- 
creases clockwise, from (cos 2 fi = 1 for perfect align- 
ment to (cos 2 6 xa } s — 0.88 for naphthalene molecules at 
1 K probed by a 100-ps x-ray pulse. By increasing the 
x-ray pulse width from 1 to 100 ps, the effective align- 
ment at T — 10 mK decreases from (cos 2 xa ) s = 0.967 

to (cos 2 9 xa ) eS = 0.944 [cf. Fig. gb) and (d)], and is 
significantly smaller than the impact of temperature rise 
from 10 mK to 1 K [cf. Fig. gc) and (d)]. 
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FIG. 6: (Color online) X-ray diffraction signal S(Q y = 0) of 
naphthalene, (a) Perfect alignment, (b) T = 10 mK, tx = 
1 ps, x-ray pulse centered at t = 20.75 ps (cf. Fig. [3|. (c) 
T — 1 K, t x — 100 ps, x-ray pulse centered at t — ps. (d) 
T = 10 mK, tx = 100 ps, x-ray pulse centered at t = ps. 
For the alignment pulses in (b), (c), and (d), we assumed 
I = 5 TW/cm 2 and t l = 100 ps. 



Although the basic features of the x-ray scattering pat- 
tern for perfectly aligned naphthalene (Fig. [6^) are pre- 
served in the scattering patterns for laser-aligned naph- 
thalene [Figs. [6^b), (c), and (d)], the contours are washed 
out and the contrast between maxima and minima is less 
pronounced with decreasing alignment. Incoherent aver- 
aging for laser-aligned naphthalene renders the diffrac- 
tion patterns more cylindrically symmetric with respect 
to the y direction, which limits the accessible structural 
information particularly at high momentum transfer. 

Side views of naphthalene in Q space are shown in Fig. 
[7| where Fig. [7^ a) depicts the long side (Q x = 0) and Fig. 



[7](b) the short side (Q z = 0) of naphthalene. Both dis- 
play strong similarities to multi-slit diffraction patterns, 
consistent with the planar structure of naphthalene being 
well aligned in the xz plane. 
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FIG. 7: (Color online) (a) Diffraction signal of naphthalene 
for Q x = 0. (b) Diffraction signal of naphthalene for Q z — 0. 
Other parameters are the same as for Fig. |6[ 



The achievable real-space resolution depends on two 
aspects. On the one hand, the smallest resolvable di- 
mension ci co h is a function of the degree of alignment [Eq. 
(43)]. On the other hand, the pixel size of the real-space 
structure reconstruction is determined by the maximum 
momentum transfer Q max in the diffraction signal S(Q). 
Here, Q max = 2tt A" 1 , corresponding to a pixel size of 
0.5 A. The approximate range in momentum space within 
which the assumption of coherent scattering holds may 
be defined by (X oh = 2n/d co h- In Fig. [8j the diffraction 
signals of Fig. |6]are shown up to Q max = 8-7T A -1 . The 
respective ranges defined by Q co h are highlighted. It is 
not possible to increase the resolution of the real-space 
structure reconstruction by choosing <5 ma x much greater 
than Q co h- In fact, numerical tests have indicated that 
data beyond Qcoh can lead to poor convergent structures. 

For 12-keV x-ray photons, a real-space pixel size of 
0.5 A requires the detection of photons scattered up to 
60° with respect to the x-ray propagation axis. The 
azimuthally averaged differential x-ray scattering cross 
section per naphthalene molecule for Q — 2n A -1 is 
da/dQ^vg = 0.6 barn. Since the largest distance between 
carbon atoms in naphthalene is ~ 5 A, the area of a pixel 
in momentum space must not exceed 47r 2 /25 A -2 . At a 
scattering angle of 60° , this corresponds to a solid angle 
c?Sl = 0.012. For a molecular beam width of 1 mm and 
an x-ray focus area of 100 /im 2 , it has been demonstrated 
that the number of molecules in the interaction volume 
can be as large as 10 7 [S5]. Hence, requiring a minimum 



of five scattered photons per pixel (115] , the estimated ac- 
quisition time for one diffraction pattern is around 70 sec- 
onds with an x-ray fluence of 10 13 photons/pulse/mm 2 
and a repetition rate of 1 kHz. 
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FIG. 8: (Color online) Same as for Fig. [6] except Q max = 
8ir A~ . The momentum transfer Q co h (see the text) is high- 
lighted by a solid green circle in (b)-(d). 
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FIG. 9: (Color online) Structure of naphthalene, recon- 
structed from the diffraction signals shown in Fig. [6] The 
pixel size is 0.5 A. (a) Perfect alignment, (b) T = 10 mK, 
tx = 1 ps, x-ray pulse centered at t — 20.75 ps. (c) T = 1 K, 
TX = 100 ps, x-ray pulse centered at t = ps. (d) T = 10 mK, 
tx = 100 ps, x-ray pulse centered at t = ps. 



C. Structure reconstruction 

We applied the phase-retrieval algorithm described in 
Sec. Ill B| to the naphthalene diffraction signals shown in 
Fig. 6 The maximum momentum transfer of 2n A -1 in 
Fig. 6 corresponds to a pixel size of 0.5 A, which is of the 
order of the structure we want to resolve in naphthalene. 
As a consequence, each pixel encodes a lot of structure 
information. Calculations have shown that the recon- 
struction routine becomes sensitive to the set of initial 
random phases. Therefore, we follow the spirit of Ref. 
[101] and average over reconstructions obtained for 100 
different sets of initial phases to define a quality criterion, 
which we apply subsequently to all 100 reconstructions to 
select the most meaningful results and average over these 
selected reconstructions. Our quality criterion is defined 
as follows: (1) All pixels of the unfiltered averaged result 
that have at least 20% of the maximum electron density 
are selected to build a density core region. (2) The qual- 
ity criterion uses the density core region and selects only 
the reconstructions that have at least 50% of their total 
electron density within this density core region. 

The filtered averages for the x-ray scattering patterns 
of Fig. [6] are displayed in Fig. [9] In our calculations, we 
employed a grid spacing in Q space of 0.16 A" 1 , corre- 
sponding to a maximum object size 8 times larger than 
the size of naphthalene. 

Figure [9] illustrates that with better effective alignment 
in the diffraction pattern more structural information can 
be reconstructed. In the perfect alignment case (cf. Fig. 
[9]ja), the positions of the atoms can be resolved. A com- 
parison with the perfect alignment reconstruction in Fig. 
|l|b), where the pixel size is 0.125 A and no averaging 
over several initial random phases was performed, illus- 



trates that the reconstruction of the carbon positions is 
more difficult when the pixel size is comparable with the 
size of the atoms. Incoherent averaging for imperfectly 
aligned molecules affects the effective resolution further, 
since the assumption of coherent scattering limits the us- 
able diffraction signal to Q < Q co h- As a consequence, 
Qmax should be of the order of Q co h and, therefore, the 
pixel size in the reconstruction is limited by the degree 
of effective alignment. 

All results displayed in Fig. [9] recover the two-ring 
structure of naphthalene, indicating an effective resolu- 
tion below 2 A. Especially the two carbon atoms con- 
necting both carbon rings are distinctly visible. Figure 
|9jb) and (9^d) show the reconstruction for naphthalene 
at 10 mK with tx = 1 ps and tx — 100 ps, respectively. 
Both have a similar effective alignment and reconstruc- 
tion. In both cases, the electron density peaks on the two- 
ring structure indicate the positions of all carbon atoms 
and are consistent with the theoretically expected effec- 
tive resolutions, which are below the typical distance be- 
tween neighboring carbon atoms (1.4 A). The predicted 
effective resolution at 1 K (c? co h = 1.73 A) is sufficient to 
resolve the overall structure of naphthalene, which can be 
seen in Fig. |9jc). For the central region of naphthalene, 
the effective resolution improves and makes it possible to 
identify the two central carbon atoms at 1 K. 



V. CONCLUSION 

We have theoretically studied the theory of x-ray 
diffraction from asymmetric-top molecules that are three- 
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dimensionally aligned by elliptically polarized light and 
have discussed a phase-retrieval method in order to re- 
construct molecular structure from the x-ray scattering 
pattern. The interaction of the gas phase molecules with 
the laser and x-ray fields was studied in the density ma- 
trix formalism. We have assumed rigid rotor molecules. 
Symmetries in the time-dependent Hamiltonian and in 
its solutions as well as in the angular density distribu- 
tion were exploited and used to improve numerical effi- 
ciency. In combination with further computational as- 
pects, a significant enhancement in numerical efficiency 
has been achieved. 

A feasibility study of x-ray diffraction from three- 
dimensionally laser-aligned molecules has been per- 
formed using the organic molecule naphthalene. We 
have linked the degree of alignment to an effective res- 
olution and have studied the impact of x-ray pulse du- 
ration and temperature on the diffraction patterns. A 
phase-retrieval reconstruction was performed on diffrac- 
tion patterns taken for different effective alignments. The 
enhancement of incoherent averaging due to imperfect 
alignment destroys structural information and hinders 
the recovery of detailed atomic configurations within the 
molecule. The reconstruction of naphthalene confirms 
our discussion that the degree of alignment is a good 
indicator for the achievable resolution in the reconstruc- 
tion. To decode structures on an atomic length scale, 
high degree of alignment in all three molecular axes must 
be achieved. Molecules have to be cooled down to a few 
Kelvin or even sub-Kelvin temperatures so rotational mo- 
tion is sufficiently reduced. Experiments have shown that 
it is feasible to cool molecules to 1 K [72"] EH [TT7j. To 
accomplish rotational temperatures well below 1 K is a 
great challenge and would open the opportunity to im- 
age large gas-phase molecules at atomic resolution. The 
problem of incoherent averaging over a finite range of dif- 
ferent orientations has to be further addressed, especially 
in the case where molecular symmetries do not coincide 
with the symmetries of the diffraction pattern imposed 
by alignment. 

With the ability to detect molecular structure, x-ray 
scattering from gas phase molecules can be used to study 
torsion effects and laser-induced deformations that are 
expected to occur during alignment in the presence of 
intense laser pulses |118j . By systematically varying the 
delay time between pump (laser) and probe (x-ray) pulse, 
it is possible to follow molecular motion on an ultrafast 
time scale. Of particular interest is the study of physical 
and chemical processes in the presence of intense laser 
fields, which simultaneously provide the required align- 
ment for imaging of reactions with atomic resolution in 
space and time. 
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Appendix A: Nuclear spin statistical weights 

In the density matrix, statistical weights, wj T = 
gj T cxp *- z(T) kT ^ i define the relative number of particles 
in an ensemble that are in a given quantum mechani- 
cal state for a well-defined temperature T. The nuclear 
spin statistical weights gj T are based on symmetry argu- 
ments and represent the number of the allowed nuclear 
spin states for a given rotational state | JtM) . 

The overall symmetry r to t of the total molecular wave 
function is independent of rotational or nuclear spin 
states and determined by the spin statistic theorem |82) . 
However, r tot is also a direct product of symmetries of 
the different quantum states [77] . 



itot — r,. 



(Al) 



where T r is the rotational state symmetry and r ns is 
the nuclear spin state symmetry. Molecules are in their 
electronic and vibrational ground states and their contri- 
butions can be omitted in Eq. ( Al I. 



The representation of symmetries depends on the sym- 
metry group of the molecule. The molecular symmetry 
group of asymmetric-top molecules without an inversion 
center is Dih or lower [119 . The rotational symmetry 
group of an asymmetric-top rotor is always V (isomor- 
phic to D2h) |1 19H121) . If there is no common symme- 
try class between r r and r ns , gj T = 1. In other words, 
molecular symmetries come only into play when molecu- 
lar rotations are identical to particle exchanges. 

All symmetry representations are written as linear 
combinations of irreducible representations (irrep) of V 
in Table Mil In addition the character table for the nu- 
clear spin and rotational states are shown. The four 
molecular rotations and particle exchanges defining V 
are: E, C a , C b , C c . 



TABLE III: Character table and irreducible representations 
of overall, rotational and nuclear spin states for naphthalene. 
The symmetry point group is V . 







Operations P of V 




irrep. 


E 


c c 


Cb Ca 


J even 


J odd 


Itot 


1 


1 


1 1 




A 


E+ 


1 


1 


(-i) j {-iy 


A 


B c 


E- 


1 


1 


(-i) J+1 M) J+1 


Be 


A 


o- 


1 


-1 


i-iy (~i) j+i 


B b 


B a 


o+ 


1 


-1 


(-i) J+1 (-i) J 


B a 


B b 


Tns 


2 8 


2 4 


2 4 2 4 


76A + 60{B a + B b + B c ) 
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The characters of the overall wave function x t ot can be 
derived from the Pauli principle j83j and correspond to 
the overall sign changes induced by the molecular rota- 
tions. First, we assume all carbon atoms in naphthalene 
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are C with nuclear spin 0. By remembering only odd 
permutations of half-integer particles (fermions) change 
the sign of the overall wave function [77], we find all 
group operation leave the overall sign for naphthalene 
unchanged, since all Hydrogen permutations are even. 
Hence, Xtot [P] = 1, VP € V and r tot = A. 

The rotational symmetry classes, in which each rota- 
tional state can be classified (cf. Sec. II A), coincide with 



the irreps of V |74j . However, whether or not J is even or 
odd defines which irrep corresponds to which symmetry 
class (cf. Table |Hl|). 

The sum of characters over all possible nuclear spin 
states Xns [P] for each operation P can be derived by the 
formula [E2J [[20] 



Xns[P] 



n (^G; 

Gp 



1) sgn(P G J 2S ^, (A2) 



where the atom group Gp is defined by the minimum 
number of identical atoms in the molecule such that af- 
ter the atom permutation P each atom in Gp is replaced 
by an atom of Gp. sgn(P<3 P ) is the sign of the permu- 
tation P within in the atom group Gp, which is +1 or 
— 1 depending whether the permutation is even or odd in 
Gp. The smallest atom group consist of one atom, which 
does not change position under the operation P. Sc P is 
the nuclear spin of the atoms in Gp. 

Imposing Eq. ( |A1| ) , nuclear spin and rotational sym- 
metries are linked together. (For naphthalene: r to t = 
T r ® r ns = A). The multiplication of symmetry repre- 
sentations is done by multiplying their characters. The 



nuclear spin statistical weights gj T are the number of ir- 
reps in r ns such that the product with a given rotational 
class is r tot . As an example: When | JtM) € 0~ and J 
is even, then T r £ Pf, and only the product Pf, <£> B\, = A 
[120] fulfills Eq. ( Al ). Since there are 60 B>b irreps in r ns , 
the nuclear spin statistical weight is 60. For all rotational 
states of naphthalene, the nuclear spin statistical weights 
are 



9Jr 



76 , r r e A 

60 , otherwise. 



(A3) 



Appendix B: Matrix elements of cos 6i m 

The squares of the matrix elements of the rotation ma- 
trix R(4>,9,x), i-e. cos 2 9i m , give a complete picture of 
the three-dimensional alignment. Knowing four cos 2 #; m 
is sufficient to describe all, where at least two axes of each 
reference frame have to be involved in the four cos 2 9i m . 
Our choice is / € {x, z},m £ {a, c}, since it is one of the 
least computationally expensive choices. The cos 2 #; m - 
matrices are diagonal in M or K when I = z or m = c, 
respectively. Only the matrix cos 2 9 xa is nondiagonal 
in K and M. Each cos 2 9i m can be written in terms 
of D^ K (cj),9,x) EH], where the matrix elements of the 
Wigner D-functions are given in Refs. [121 IMl HI] ■ The 
matrix elements of the four matrices that we have chosen 
are: 



(JKM\ cos 2 9 ZC \J'K'M') 
(JKM\ cos 2 9 XC \J'K'M'} 



(JKM\ cos 2 9 za \J'K'M'^ 



(JKM\ cos 2 9 xa \J'K'M') = 



l s , . 2 2J+1 

3 6jj,6kk ' 5mm ' + 2FTi 

x ( J,M; 2, 0| J',M') (J,K;2,0\J',K') , 
I . . II / 2.7 - 1 



;$JJ'$KK<SmM' 



6 V 2 J' + 1 



J,A";2,0| J',K' 



(Bla) 
(Bib) 




! \l -(J,M; 2, 0| J',M') - ( J, M; 2, 2| J , M ) + ( J, M; 2, —2\ J , M ) 



1 /2J + 1 



6 V 2 J' + 1 




; JM; 20 1 J'M') 



(Blc) 



J,K;2,2\J',K') + (J,K;2,-2\J',K') 



l s , . 1 2J+1 

T jj ' 6kk ' Smm ' + IV 2STl 



< I W - { J, M; 2, 0| J , M ) - ( J, M; 2, 2\ J , M ) + ( J, M; 2, -2| J , M 



(Bid) 




M t/g<J s A';2,0|J / s A")- 



;j,K;2,2\ J>, K') + (J,K; 2,-2| J',K' 
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